Local entropic effects of polymers grafted to soft interfaces. 
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In this paper, we study the equilibrium properties 
of polymer chains end-tethered to a fluid membrane. 
The loss of conformational entropy of the polymer 
results in an inhomogeneous pressure field that we 
calculate for gaussian chains. We estimate the ef- 
fects of excluded volume through a relation between 
pressure and concentration. Under the polymer pres- 
sure, a soft surface will deform. We calculate the 
deformation profile for a fluid membrane and show 
that close to the grafting point, this profile assumes a 
cone-like shape, independently of the boundary con- 
ditions. Interactions between different polymers are 
also mediated by the membrane deformation. This 
pair-additive potential is attractive for chains grafted 
on the same side of the membrane and repulsive oth- 
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I. INTRODUCTION. 

Fluid membranes are surfactant bilayers self-assembled 
from solution jjj . They are the prevalent constituents of 
many natural and industrial colloidal suspensions, that 
often contain also other macromolecular species. In the 
biological realm, phospholipid bilayers build the walls of 
liposomes and cells, hosting proteins responsible for func- 
tions as diverse as anchoring the cytoskeleton, providing 
coating protection against the body immune response or 
opening ionic channels for osmotic compensation j^j. In 
cosmetics, pharmaceutics or detergency, many formula- 
tions are membrane solutions with polymers added for 
performance, processing, conditioning or delivery ||. 

The interactions between polymers and fluid bilayers 
have been well scrutinized in many systems. Polymers 
grafted to the bilayers can induce gelation Q or other 
phase changes in liquid lamellar phases. They sta- 
bilize monodisperse vesicles and modify the geome- 
try of monolamellar (8||] and multilamellar cylindrical 
vesicles |l(J. Theoretically, the behaviour of fluid mem- 



branes is well understood in terms of bending elastic- 
ity pyr^ ], a description that requires as an input the 
value of three material constants: the bending rigidity 
k, the gaussian rigidity R and the spontaneous curvature 
radius Rq. One might hope that the behaviour of mixed 
systems can still be described by an effective elastic en- 
ergy, with modified material constants. The task that 
theoretical studies have undertaken is to calculate the 
modifications induced on k, R and Rq by the addition of 
the macromolecules p^-[l7|l . 

However, polymer-membrane interactions must have a 
local quality. For instance, if a polymer is end-tethered 
to a membrane, it is clear that the interactions are strong 
at the anchoring point and vanish far enough from it. We 
show in this paper that, for grafted polymers, it is possi- 
ble to construct a local description of polymer membrane 
interactions. Our description stands on the recognition 
that an end-grafted polymer applies a pressure field to 
the grafting wall 18 llj. The pressure field results into 
a local deformation of the membrane: a membrane with 
grafted polymers can therefore be seen as a surface with 
bending elasticity carrying a number of pressure patches, 
each of them creating its own deformation. 

The paper is organized as follows. In the next section 
we compute the pressure field applied by a grafted poly- 



mer in theta and good solvent conditions. In section III 
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we consider the case of a freely-standing membrane, for 
which we compute the deformation induced by the poly- 
mer pressure patch. We also show in this section that 
the interactions between the different deformation fields 
give rise to a membrane mediated potential between dif- 
ferent grafted polymers. The consequences of such po- 
tential are briefly discussed. Section IV is dedicated to 
two membrane geometries relevant for experiments. We 
first discuss the case of supported membranes, and then 
the case of lamellar phases. In the conclusions, we will 
briefly speculate on the implications of our results for 
hairy vesicles. 



II. THE PRESSURE APPLIED BY A GRAFTED 
POLYMER. 

The number of available conformations for a long, flex- 
ible polymer is strongly reduced by the process of end- 
grafting the chain to a hard wall 20 , 2l[[. The average 
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configuration of the macromolecule is a compromise be- 
tween the need to avoid the surface and the constraint 
imposed by the tethered end. It is clear that if the sur- 
face can be deformed, there will be an entropic reason 
to push it away from the monomer cloud. This can be 
described as a pressure that the polymer applies to the 
grafting wall. In the following paragraph we explicitly 
compute the pressure for ideal chains and relate it to the 
concentration at the wall. We then argue that this also 
provides a good pathway to compute the pressure applied 
by chains with excluded volume. 



A. Gaussian grafted chain. 

We consider a gaussian chain of N units, end-tethered 
by one extremity to a non-adsorbing wall. The surface is 
described by its height h(x,y), where (x,y) denotes the 
position in the horizontal coordinates frame. The ther- 
modynamic properties of the chain are described by the 
propagator Gn that satisfies the Edwards equa- 

tion psji 



dG N (f,f') a 2 
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AGjv(r,f) 



(1) 



with the following boundary conditions: Gn{^,t') = 
on the wall and limjy^ oG\zv(r, r') = 5(f — r'). The 
length a is the monomer size. The statistical weigth of 
the chain attached at a monomer distance from the origin 
a = (0, 0, a) is given by 



Z N (a) = / df'G N (a,f') 



(2) 



the integral running over all the space available for the 
free end. In the flat, reference case h(x, y) = 0, the Green 
function can be factorized 



term of the perturbative expansion obeys an Edwards 
equation 



dZ 



(0 2 
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The solutions of successive orders are coupled through 
the boundary conditions on the wall 



= Z N (x,y, h) 

t(x,V,Q) 

h 2 (x,y) d 2 Z N 
2 dz 2 



d Z 

Z N (x,y,0) + h(x,y)-^—(x,y,Q) 



(x,y,0) + 
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In the following, we concentrate on the first order term 
ZjPi that is related, as we will see, to the pressure ex- 
erted by the polymer on the surface. Zff is given by the 
solution of equation (||) with the boundary condition 

Z$\x,y,0) = -h(x,y)?^-(x,y,a) (7) 



The solution can then be written as []2Jj 

-f(O) 
r _N_ 

Or: 
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so that the displacement of the surface is achieved to first 
order in h by the work 



AF = F[h] - F[0] 
= -fc fl Tlog[l 



^jv - 

7 (0) J 



= / dSp(x,y)h(x,y) 
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and the partition function is 
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where the function p(x, y) has the radially symetric form 
k B T 
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with R g — yjNa 2 /6 the gyration radius of the chain, and 
erf the error function pjjfl . Now we seek for a perturba- 
tive solution p3] of the Edwards equation by performing 
a small displacement h of the surface. We write the parti- 



tion function as Z 



N 



is of order h l and zffi is defined in equation (k 
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where Z 
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can notice that due to the linearity of equation (111), each 



with r = y 1 x 2 + y 2 . At point f = (x,y), the elementary 
work dF required to displace a volume dV(r) = h(r)dS is 
given by dF — p(r)h(r)dS . The function p(r) is therefore 
the pressure applied by the polymer on the surface at 
point r. It is a non-homogeneous function - see figure (|l|) 
- that sharply decays from its maximum value at the 
anchoring point with the scaling form 



p(r) ~ - — — for a < r < R„ 
27rr J 



(11) 



Well inside the polymer umbrella (a <C r <C R g ), where 
the only relevant length is r, expression ( |TT| ) is the nat- 
ural scaling for the pressure. In this region, most of the 
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monomers that contribute to the pressure are close to the 
grafted end. For distances larger than the polymer size, 
the pressure vanishes exponentially. A grafted polymer 
can then be pictured as a microscopic pressure tool that 
applies a well-defined but non-homogeneous force on a 
disk of radius ~ 2R g centered at the anchoring position. 
In the middle of the patch the pressure has a strong value: 
p(0) ~ 2.4 x 10 7 Pa for a = 0.3nm at room temperature 
T = 25°C. A small area within a monomer distance from 
the origin supports most of the total force / exerted by 
the chain onto the surface 



1 h T ^ 

2irrdrp(r) = -Jl_exp(-^) ~ 13.3pN (12) 



with the precedent values of monomer size and tempera- 
ture. The grafted monomer exerts a point like force — / 
that ensures mechanical equilibrium. 

Previous work (ll]]2(| has focused on curvature con- 
tributions to the polymer free energy AF: by consid- 
ering a surface of a given shape, AF is calculated as 
a function of the curvature Jj. For instance, for a 

sphere and a cylinder one gets AF sp h — —^UbT^ 
and AF cy i = ^AF sp h- The minus sign indicates that 
the energy is lowered by spherical and cylindrical sur- 
faces that bend away from the polymer. We naturally 
recover these results by considering a general surface de- 
fined by h(x,y) = — — 2B^' an( ^ evaluating the inte- 
gral (H) with Ri = R 2 = R for a sphere and R± = R , 
i?2 = for a cylinder. Interestingly, for a minimal sur- 
face = — i?2 = R), there is no contribution of the 
curvature to the free energy (AF — 0): to first order, the 
cntropic cost of tethering a gaussian chain to a plane or 
to a minimal surface is the same. 



B. Relation between pressure and concentration. 

As explained above, the pressure that a grafted poly- 
mer exerts on the wall has an entropic origin: by displac- 
ing the surface at point f = (x,y) from its flat position 
h(f) one increases (h < 0) or decreases (h > 0) the num- 
ber of allowed chain configurations. The work per unit 
surface associated with the corresponding entropy gain 
or loss defines the pressure. Alternatively, the pressure 
can be viewed as resulting from the forces applied by all 
the monomers at a given surface point. If the surface 
potential acting on each monomer is u(z), the pressure 
is given by [p?! 



p(x,y) = - 



dz^-(z)c(x,y,z) 



(13) 



with c the monomer concentration. Expression ( |13| ) re- 
veals the linear relationship between pressure and con- 
centration at the wall but is not very usefull for the con- 
tinuous gaussian chain considered in this paper. Instead, 



we directly derive the pressure-concentration relation by 
noting that the monomer concentration of a chain teth- 
ered to a flat surface is written as M 
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The interaction with the wall being purely repulsive, one 
has c^°^(x,y,0) = and ^^-(x,y, 0) = 0, the second 
derivative of the concentration being the lowest derivative 
that does not vanish on the surface. Rewriting equation 
(@) by taking into account the definition of the partition 
function (0) leads to 



a 2 9 2 c (0) 



(15) 



with r — \J x 1 + y 2 . Qualitatively, the pressure can be 
associated with an ideal gaz pressure caused by the con- 
centration of monomers at a distance z = -4= from the 

\/6 

wall 



p(r) = k B TcW(r,z = — ) 



(16) 



Equivalently, equation ( |13[ ) can be used to assert that the 
effective wall potential acting on the monomers has a sec- 
ond moment of forces given by / °° dzz 2 ^(z) = ksT^-. 



C. Grafted self-avoinding walks. 

Flexible polymer chains in good solvent cannot be de- 
scribed by gaussian statistics. They only exhibit ideal 
gaussian behaviour close to the theta temperature, at 
which monomer attraction compensates for steric repul- 
sion. Above the theta point, the polymers perform self- 
avoiding walks (SAW's) which lead to distinct statistics. 
In particular the average dimension of a SAW coil is 
larger than its gaussian equivalent, the end-to-end dis- 
tance scales as R — N"a, with v the Flory exponent, 
close to v ~ 3/5. 

As stated in equation ([l3]) , the proportionality between 
the polymer pressure and the monomer concentration in 
the vicinty of the wall is expected to hold on general 
grounds, independently of the approximations involved. 
Altough equation (|l^) does not easily provide for a pro- 
portionality coefficient, it does give a firm ground for 
predicting the scaling form of the pressure applied to 
the wall by grafted chains in a good solvent. For com- 
paraison, we first recall the structure of the monomer 
concentration profile for a gaussian grafted polymer, a 
case where an explicit calculation can be performed | pl| . 
The cone of equation z = r separates two regions in 
space. Outside the cone (x<r) but well inside the poly- 
mer "umbrella" (z, r <C Rg), the monomer concentration 
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grows quadratically from the wall, c(r, z) ~ z 2 /(r 3 a 2 ). 
In this region the wall has an important depletion effect 
on the polymer configurations. Inside the cone (z 3> r) 
one recovers the usual bulk concentration of a gaus- 
sian chain, c(r) ~ \/(ra 2 ). The crossover between the 
two behaviours occurs as one crosses the cone surface 
z = r. In the presence of excluded volume interactions, 
the bulk concentration is given by c(r) ~ l/(r 4 / 3 a 5 / 3 ). 
For chains with excluded volume, the profile grows from 
the wall as c(r,z) ~ z 5 / 3 . Writing the scaling form 
in the region z <C r that matches bulk behaviour at 
z = r leads to c(r,z) ~ z 5 / 3 /(r 3 a 5 / 3 ). The pressure 
applied by a swollen grafted chain is therefore given by 
p(r) ~ fcflTc(r, z = a) ~ fc s Tr~ 3 . It has the same scaling 
form as the pressure applied by ideal chains. We expect 
such scaling to be rather independent from the molecular 
details or from the differences between chain models. 

Excluded volume effects might nevertheless modify the 
amplitude and the range of the applied forces. In order 
to quantify these effects we implemented a Monte-Carlo 
simulation on a polymer attached to a flat, impenetrable 
wall. The chain is described as a pearl necklace ^8|, each 
pearl of size a. The first monomer is grafted to the wall 
with center of mass position (0, 0, 0). During simulation a 
histogram c(r, z) for the monomer center of mass concen- 
tration at a distance r along the wall and height z above 
the wall is compiled. We are using a binsize of 0.18 • a 
in the r and z direction. To extract the concentration at 
the wall at a distance r away from the grafting point we 
fit the function f(z) = c(r, z) with a fourth order poly- 
nomial in z multiplied with an exponential exp(— Xz), X 
a fitting parameter. Upon extrapolating the fitted func- 
tion to z = we get the concentration of monomers at 
the wall at a given distance r: lim z _»o c ( r : z )- To ensure 
a reasonable error bar on the resulting concentration we 
generated 6T0 6 configurations for gaussian random walks 
and self-avoiding random walks. By analysing the statis- 
tics for end-to-end distance, which represents the slowest 
relaxing mode for the chain, we estimate that we have an 
maximal error-bar of 12% as r is increased from a few a, 
where the error-bar is less, to 30 • a. 

We first consider the Monte-Carlo results for a chain 
without excluded volume. Because our Monte-Carlo 
chain is actually a freely-hinged chain, we do not expect 
the amplitude coefficient of the gaussian model to ex- 
actly hold. We therefore plot in figure (0a), both the ex- 
pression for the pressure from equation (|10|) and the val- 
ues for the monomer concentration at the wall extracted 
from a Monte-Carlo simulation of a chain with N = 200 
monomers. Agreement is excellent, except at distances 
of the order of a monomer size where the fixed length be- 
tween monomers induces oscillations reminiscent of those 
observed in the correlation function of hard spheres. In 
figure (^b), we show equivalent results for a chain of 
N= 200 monomers with excluded volume. Agreement 
is also excellent, if we replace the dimension of the chain 



R 2 = No 2 in expression @ by R 2 = l.5N 2u a?, with 
v = 0.6. For the chain representation used in our simu- 
lations, excluded volume effects influence only the range 
of the pressure field. The scaling form at small distances 
and even the amplitudes are equivalent to those of the 
ideal chains. 

III. POLYMERS ANCHORED ON A FREELY 
STANDING MEMBRANE. 

Grafted polymers are small pressure patches: when a 
polymer is grafted to a soft interface, the pressure de- 
forms the interface into a characteristic shape which is 
determined by the balance between the pressure and the 
elastic response of the grafting surface. In this chapter we 
consider first the deformation induced by a chain grafted 
on a freely standing elastic membrane, and then mem- 
brane induced interactions between two grafted chains. 

A. Deformation induced by the pressure field of a 
gaussian polymer. 

The thermodynamic properties of a fluid membrane 
are well described by the Canham-Hclfrich Hamilto- 
nian |5L~3j| a provided that the thickness d of the bilayer 
is small compared to the other relevant lengths of the 
problem (i.e. d<^R g ). The surface is described by its 
height h(x,y), where r= (x,y) refers to the coordinate 
frame in the reference plane h(x,y) — 0. Assuming a 
gentle surface deformation (| V/i |<C 1), the Hamiltonian 
is written in the Monge representation as 

E =\\ dxdy{Ah) 2 (17) 

with k the bending rigidity and A the 2-dimensional 
Laplacian operator. For a membrane with fixed topol- 
ogy, the Gauss-Bonnet theorem implies that the gaussian 
curvature term is constant and can therefore be ignored. 
For a purely repulsive surface, the total free energy of the 
system {membrane-!- polymer} is the sum of the bending 
energy of the membrane and the work of the entropic 
force exerted by the polymer 

F[h] = F[0] + ~ J 2nrdr(A r h) 2 + J 2Ttrdrp{r)h{r) 

(18) 

with F[0] the work required to graft the polymer to a flat 
plane. Since both the applied pressure and the boundary 
conditions considered below are radially symmetric, we 
coveniently expressed all quantities in cylindrical coordi- 
nates. The equilibrium shape of the membrane ensues 
from a compromise between the applied pressure and the 
restoring bending forces. Functional minimization of the 
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free energy with respect to the membrane profile h(r) 
leads to the Euler-Lagrange equation for the equilibrium 
shape 



K,A r A r h{r) +p(r) = 



(19) 



where A r = ±£r£ 



is the radial part of the Laplacian 
operator. We first focus on the central region, close to the 
grafting point, where most of the stress is concentrated. 
Here (r <C R g ), the pressure behaves likes /csT/(27rr 3 ), 
so that the resulting deformation has a cone shape 



h(r) 



, k B T 



2tt 



(20) 



independently of the boundary conditions. A point- 
like defect is generated at the origin and we will refer 
to this conic shape as the fundamental pinch - see fig- 
ure ([?]). The surface curvature diverges at short dis- 
tances as Arhtxr^ 1 . Physically, this divergence is cut- 
off either at a distance of the order of the membrane 
thickness d by non-harmonic terms in the curvature en- 
ergy or, for infinitely thin membranes, at the monomer 
length a at which the pressure saturates. It is inter- 
esting to note that if a polymer is grafted to the tip 
of a purely conic deformation with an arbritary slope, 
and then the slope determined by balancing the global 
entropy gain of the polymer and the elastic cost of de- 
forming the membrane into a cone shape, one finds the 
same slope as that of expression (p0|), given by the lo- 
cal balance of equation ([!{]) Jl(|. This is an indica- 
tion that the conic region supports most of the total 
stress imposed by the pressure patch. Neglecting ef- 
fects in the cut-off region, the analytical solution of equa- 
tion ( |l9|) can be written as h(r) = h p (r) + hbh(r), with 
hbh(r) — [c\ + c 2 ln(r) + c 3 r 2 + c 4 r 2 ln(r)] the kernel 
of the biharmonic operator and 
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'2^ [ r exp( -4i?f' 



R, 



■erfc( 



2R 



4 9 y 2R g ' 



Rn 



* -<)] (21) 
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a particuliar solution of equation ( |l9[ ) . The constants C\ , 
c 2 , C3 and C4 are determined by the boundary conditions. 
In the simple case considered here, the membrane has 
no imposed constraints other than its known position of 
the center of coordinates h(0) = 0. Because there are no 
forces acting on the membrane at large distances from the 
center the average curvature must vanish there A r h(r — > 
00) = 0. This determines the four constants c\ = C2 = 
C3 = C4 = 0. At distances larger than R g , the profile is 
a catenoid, a radially symetric shape with zero average 
curvature 



MO 



k B T 
2ttk 



R9H-5-) fovr^R g 
K 



(22) 



The complete profile is displayed in figure (§. It has 
a characteristic pinched shape, with a cone like defor- 
mation (|l^) that crosses over to the catenoidal shape 
( p2| ) . The divergence of the profile is related to the un- 
constrained nature of the membrane considered here. We 



will discuss in section IV how the deformation profile is 



modified by the boundary conditions or other external 
fields. 



B. Interaction potential between two grafted 
polymers. 

Most often, bilayers anchor a finite concentration of 
polymers. Each polymer is a pressure patch that car- 
ries with it a pinched form. Beyond the usual Van der 
Waals or steric interactions between the different chains, 
the superimposition of the different pinches will also lead 
to membrane mediated forces. Due to the linear nature 
of the pressure contribution to the free-energy of the sys- 
tem, the many-body problem reduces in this case to a 
sum of two body interactions that we now study. 

The free energy of two pressure patches applied at po- 
sitions r± and r*2, on the same side of a membrane is 

F[h,n,f 2 ] = I dS\p{\ r-n \)+p(\ r-r 2 \)]h{r) 



dS(Ah(r)) 2 



so that the deformation field obeys 



(23) 



KAAh(r)+p(\ f-ri |) +p(\ r-r 2 \) = (24) 

with A = A r + -pz-§ffi the Laplacian operator in cylin- 
drical coordinates. Note however that the problem has 
now lost its radial symmetry. A particular solution of 
this equation is h p (\ f—fi |) + h p (\ r — r 2 |), the func- 
tion hp being given by (|2l]). The general solution of the 
biharmonic equation AAh = that satisfies the require- 
ment lim r ^ oo Ah = with a finite value at the ori- 
gin is simply a constant. If we impose the conditions 
h(fi) — h(r 2 ) — 0, we are lead to the solution of the 
differential equation ( |2^ ) 



h{r) = h p (\ r- n |) + h p (\ r- r 2 |) - h p (l) - h p (0) 



(25) 



with I the distance between polymers, I =| f± — \. 
The interaction potential V(l) is given by the difference 
between the free energy ( p3| ) and twice the free energy of 
one isolated polymer. Inserting the solution of the Euler- 



Lagrange equation in (23) and integrating by parts leads 
to 



V curv (l) = —K I dSAh p (r)Ah p (f - I) 



(26) 
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The potential V(l) is always negative: the interaction be- 
tween two polymers attached on the same side of the bi- 
layer is attractive. It also follows from equation (|2^) that 
two polymers anchored to the opposite side of a mem- 
brane repel each other. The interaction potential has in 
this case the same functional form but with the reverse 
sign. The interaction potential is shown in figure (Q) for 
two polymers grafted on the same side. Notice that the 
potential range is of the order of the polymer size: at 
larger distances, the deformation field having a zero cur- 
vature shape, the cost of grafting a second polymer to 
it is the same as grafting a polymer to a flat interface. 
The mechanisms responsible for attraction or repulsion 
are easy to understand. When two polymers are grafted 
to the same side of a membrane, they can both share the 
same deformation profile, instead of creating each one 
a profile of their own. When they are grafted to oppo- 
site sides, the deformations are mutually neutralized: the 
polymers can only fully develop their deformation profiles 
at a large distance from each other. 

At short distances (I -C R g ), the attraction has a log- 
arithmic behaviour 

UW^MW (27) 

For soft membranes the elastic constant is of order ksT, 
which leads, for a polymer with a radius of gyration of 
ten nanometers, R g — lOnm and a minimum approaching 
distance of the order of a monomer size I = 0.3nm, to an 
attraction well of a couple of kgT. Moreover, the poten- 
tial varies quadratically with temperature, this class of 
interactions is thus quite sensitive to temperature vari- 
ations. We will further discuss in the conclusions the 
possible implications of such sensitivity. 

It is important to stress the differences between the 
potential that we just described and other membrane in- 
duced interactions abundantly described recently |^9|,^o| . 
The pressure patches do not lead to Casimir-like inter- 
actions, there are here no logharitmic or algebraic tails. 
In the literature, the inclusions are considered as rigid 
bodies that fixes the value of the membrane curvature at 
the inclusion site. It can be checked that, if one consider 
very soft point-like inclusions for which the curvature self- 
adjusts to a prefered value that minimizes the inclusion 
distortion and the membrane bending, then there is no 
interaction potential between the inclusions |^T| . In the 
jargon of membrane induced interactions, our potential 
corresponds to a "short range" potential. Notice however 
that for polymers, the range of the interaction can be at 
least one order of magnitude larger (~ tens of nanome- 
ters) than typical inclusion sizes (~ few nanometers). 

The potential ( |26| ) accounts only for curvature- 
mediated interactions. At distances of order R gi inter- 
chain interactions give also rise to a repulsive contribu- 
tion. We qualitatively account for the polymer-polymer 
repulsion by separating the two chains with a mid-plane 



hard wall J24j. This over-estimated repulsive potential 
reads — kBT\n[er{(l/R g )]. We now write the total inter- 
action potential as 

V(l) = V curv {l) -bk B Tln[etf(l/R g )] (28) 

with b a constant smaller than unity. Both parts behave 
logarithmically at short distances: the potential is attrac- 
tive for small values of the bending rigidity and is repul- 
sive for high values. Figure (|5|) shows the plots of V(l) for 
different values of the rigidity k and b arbitrarily fixed at 
b = 1/(2%). In this paticular case the crossover between 
attraction and repulsion occurs at n — kgT. Chang- 
ing the value of b will accordingly rescale the crossover 
value. When the surface is covered with a finite density 
of chains, the onset of aggregation can be monitored by 
the second virial coefficient B 

B = \J -exp{-^}) (29) 

The plot of B as a function of the bending rigidity with 
our particular choice of b is shown on figure (§: for k < 
0.6fcsT, the second virial coefficient becomes negative, 
indicating aggregation of different chains. 

C. Star-like polymer aggregates. 

Attractive interactions between chains grafted to the 
surface may lead to a star-like structure: the grafting 
points merge into a core, while excluded volume repul- 
sions acting on the arms give the aggregate the hemi- 
spherical shape displayed in figure (0). We now discuss 
how the structure of the aggregate changes the nature 
of the pressure applied to the grafting surface, and the 
deformation that the pressure induces on a free standing 
elastic membrane. 

The structure of a star polymer with / arms of poly- 
merization index N can be described by the Daoud- 
Cotton model |3^] . Attachment of the chains to a central 
core effectively forces the local polymer density to be 
everywhere inside the star above overlapping concentra- 
tion. The star can therefore be described as a semi-dilute 
solution p0| , with a local, position dependent screening 
length £(r), where r is here the distance from the center 
of the star in a frame of spherical coordinates. Pictorially, 
we represent this by associating with each arm a string 
of blobs of increasing size £(r). The radial dependence of 
the blob size £(r) can be obtained by noticing that at a 
distance r from the center there are / blobs of cross sec- 
tion £(r) 2 occupying a total area of 4irr 2 . The blob size 
thus varies as £(r) ~ rf~ 1 / 2 and the corresponding local 
polymer concentration as c s (r) ~ j 2 / 3 a - 5 / 3 r - 4 / 3 . Note 
that there is a crowded region of size R c ~ a/ 1 / 2 in the 
middle of the star where the concentration reaches one. 



G 



The size of the star can be obtained from monomer con- 
servation Nf = Air J q r 2 drc s (r) . Neglecting the small 
core region one 

getsi?~aiV 3 / 5 / 1/5 - 
The structure of interest to us is a half-star, simply 
obtained from the Daoud and Cotton model by replacing 
f/2 arms by a repulsive half plane - see figure (0). The 
local pressure is a simple function of the blob size p s (r) ~ 
ksT(,(r)~ 3 , leading to a pressure applied on the surface 
of the form 



Ps(r) c fl 2 



(30) 



if we now revert back to the previous notation where r 
is the distance from the center in a surface cylindrical 
frame. The reasons beyond this functional form can also 
be related to the interfacial structure of the star. Close to 
the surface there is a monomer depletion layer, growing 
from the surface as c(r, z) <~ z 5 / 3 . Due to screening, the 
bulk behaviour is recovered when one crosses the cone 
surface z = rf~ 1/>2 . Writting again the r dependence 
of the depletion layer such as to match the bulk value 
c s (r) ~ / 2 / 3 a ~ 5 / 3 r ~ 4 / 3 , one gets c(r,z) — p s (r)(z/a) 5 ^ 3 , 
with p s {r) the pressure field in equation (|30|). For dis- 
tances r larger than the star size, we expect the pressure 
to vanish rapidly. 

For ideal chains, the pressure applied by / chains 
grafted at the same point is / times larger than the pres- 
sure applied by a single chain. For polymers with ex- 
cluded volume there is an additional crowding effect that 
results in a pressure / 3 / 2 times larger than the pressure 
of a single chain. Also, the range of the pressure grows as 
the star size, and is a factor J 1 / 5 larger than the range of 
a single chain. Since the pressure field close to the graft- 
ing point has the same scaling form as a single chain 
pressure, the patch still induces a conic deformation on 
a free standing elastic membrane. The angle of the cone 
is more pronounced for a star than the angle of a pinch 
from a single chain 



Mr) ^ -f /2 ( — )r 

K 



(31) 



and the energy gained by the creation of the star-pinch 
is much greater than the the sum of the energy gains of 
/ individual pinches 



k H T 



k a 



(32) 



The Daoud-Cotton model also allows to compute the ex- 
cluded volume cost to build a star-like structure, the re- 
sult is 



^~/ 3/2 M-) 

Kb J- a 



(33) 



For a number of chains / greater than a threshold fo ~ 
(At/feeT) 2 / 3 , aggregation is always favored. Nevertheless, 



that process might be slow since it is hindered by an 
energy barrier AF ~ k. It is also important to stress that 
the mecanical constraints on the bilayer may limit the 
maximum aggregation number. This can be determined 
by setting | V/i |~ 1 in expression ( |3l| ) which leads to 
/max ~ fo- At this stage a piece of membrane decorated 
with / m ax polymers might as well detach from the main 
membrane, leading to a coexistance between decorated 
membranes and decorated small vesicles or micelles fij]. 



IV. POLYMERS ANCHORED ON A 
CONSTRAINED MEMBRANE. 

In most practical situations, the bilayer does not stand 
free in the solvent but it is subjected to additional con- 
straints. In this section, we consider two important situ- 
ations. First we discuss the pressure applied by a grafted 
polymer to a supported bilayer. Supported bilayers are 
fluid membranes that adhere to a substrate, allowing for 
instance for the observation of cell phenomena like mem- 
brane protein aggregation, opening of ionic channels and 
others |^a-B5j. Adhesion of membranes on a substrate 
plays a role in biological phenomena such as endocytosis 
and exocytosis Q , it is also of relevance in biotechnolog- 
ical processes, such as drug delivery by liposomes [ [36| . 
A second important case of constraint membranes cor- 
responds to lamellar L a phases, where the membranes 
are confined by interactions with their neighbors in the 
lamella stack fE|]37 



A. Supported bilayers. 

We consider a membrane adhering to a flat surface 
with a contact energy per unit area T/2. Application 
of a pressure patch peels off a region of radius L and 
central height ho = h(r = 0). The corresponding free 
energy functional is 

F[h(r),L, ho] = F - TttL 2 + - 2irrdr(A r h) 2 

2 Jq 



+ / 2nrdrp(r)(h{r) - ho) (34) 
Jo 

with Fo a constant. Notice that the membrane height h 
is here measured with respect to the substrate. Perform- 
ing the functional minimization of equation (|^) with re- 
spect to the deformation h leads to an Euler-Lagrange 
differential equation identical to equation (|nj) of the 
free-standing case, with boundary conditions h(0) = ho, 
h(L) = an d jr (L) — 0. Further minimization of the 
free energy p4) with respect to the size of the peeled 
region L, and with respect to the heigth at the origin ho, 
provides two more boundary conditions 
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dF 



D 



(-) 1/2 

Hi 



" = Q^ Ah{r = L ) = 

oho dr 



(35) 
(36) 



Relation (35) equates the balance between the attractive 



potential and the elastic moment at the contact line r = 
L |]38f . The peeled radius is given by the implicit equation 
for L 



1 



L 



1 ^P(-T^) - V^rf(— ) + *1T0 = (37) 



2 L 



2R 



Ra 



The parameter (3 2 — nTR 2 g / (hsT) 2 controls the value 
of L. For small adhesive energies or very flexible mem- 
branes we have (3 <C 1, the membrane is loosely at- 
tached to the surface and the deformation is similar to 
the free-standing case: L = i? g /(2 v /7r/3) 1 / 2 and ho = 
(kBT/iy^n^Rg \n(L/R g ). For large adhesive energies 
or stiff membranes, /3 > 1, only the conical deforma- 
tion survives and the pinch height is proportional to the 
peeled radius: L = R g /(2-K(3) and ho — (3/csT/167r«:)L. 
One migth be astonished that these quantities cannot be 
expressed only as a function of the lenght A = (ft/T) 1 / 2 . 
For many problems involving bending energies and ad- 
hesion (or interfacial tension) of membranes this length 
separates two regimes. On length scales larger than Ao, 
adhesion or tension effects dominate the behaviour of the 
membrane, while for lenghts below Ao the deformation 
is ruled by the bending curvature. In our case, the re- 
sults can also be understood in terms of Ao, by recalling 
first the simpler case of a membrane that adheres to the 
surface but has a fixed, given heigth (o at the origin. 
It is easy to show that in that case the peeling radius 
is given by L — (8A^Co) 1/4 ~ (A Co) 1/2 - But when a 
polymer patch is applied, the heigth is fixed by the poly- 
mer pressure and it follows, for strong adhesions, the 
conic form (20). One has thus Co ~ L/k which leads 
to L ~ A/k. This holds up to a length L of the or- 
der of the radius of gyration R g . For smaller adhesion 
strengths the balance is determined by Co ~ Rg/ K leading 
toL~ {RgXo/n) 1 ' 2 . 

Adhesion energies can be found in the range T ~ 
1CP 7 — 1CP 4 mN.m -1 . For typical bending modulii 
k ~ 5 — 20/cbT, the extension of the peeled zones is in the 
range of one to ten nanometers, which is also the typical 
size for polymers. 



B. Membrane in a potential well. 

In this section, we focus on the effect of pressure 
patches applied to membranes confined in an harmonic 
potential. This is for instance relevant to describe the 
lyotropic smectic phases L a but serves also as a 
paradigm for other situations where the membrane is con- 
straint by an external soft potential. The energy func- 
tional of one bilayer reads 



F[h(r)) = F[0] + / 2Trrdrp{r){h{r) - h ) 
Jo 

k f°° B f°° 

+- / 2irrdr(A r h) 2 + — / 2irrdrh(r) 2 (38) 
2 Jo 2 J 

where ho is the value of the deformation at the origin. 
In the case of a stack of membranes the amplitude of 
the harmonic potential B is the compression modulus of 
the system and can be related to the curvature of the 
interaction potential. The natural length arising in ex- 
pression dH) is = ("g) 1 ^ 4 - O n scales larger than l 
the deformation is controled by the harmonic potential 
while on smaller scales the behaviour of the membrane 
is ruled by curvature effects. In particular, for non-ionic 
L a phases where steric Helfrich interactions control the 
repulsion between the membranes |59],[l0|, it is possible 
to show that the lenght is proportional to the average 
interlamellar spacing d 



(39) 



Minimizing the free energy functional with respect to C 
h — ho leads to 



«A r A r CM + BC(r) = -(p(r) + Bh ) 



(40) 



with boundary conditions are £(0) = and C(°°) = ~ ^o- 
In order to solve equation ([l0|) ; we first evaluate the as- 
sociated Green function g(r,r') defined by 



nA r A r g(r, r') + Bg(r, r') = S(r - r') 



(41) 



so that the deformation field £ is eventually obtained 
from 



CM 



dr'g(r 7 r')(p(r)+Bho) 



(42) 



A straightforward but tedius calculation gives for r > r' 



9+{r, r') 



and for r < 



9-(r, r') 



-M^ V{Beio(ryZo)Kero(r/Z ) 

K 

+Ber (r'/7 )Kei (r/Zo) (43) 
+-Keio(r7Zo)Keio(r/Z )} 



_M^V{Bei (r/Zo)Ker (ryZo) 
+Bero(r/Zo)Keio(r'/7o) (44) 
+-Keio(r/Z )Keio(r'Ao)} 



where the Kelvin functions Ber (cc), Beio(x), Ker (x) and 
Keio (x) are the real and imaginary parts of the modified 
Bessel functions I (xe 47r / 4 ) and K (a;e I7r / 4 ) Q. The 
harmonic potential allows then for an oscillatory profile, 
as depicted on figure (1st) . 
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The amplitude of the deformation, ho is obtained by 
minimizing the free energy with respect to ho, leading to 

, k B T f°° dx, 4 , , x 1 w x 1 x 

(45) 

where a = Rg/lo ■ Asymptotically h varies like the 
logarithm of a for small values of a, and decays with a 
power law for large values: ho ~ or" 1 - see figure ([)]). For 
Helfrich systems where the length Iq is proportional to 
the interlamellar spacing Iq ~ d (and k ~ ksT), one may 
easily induce a deformation of O.ld, by using polymers 
of gyration radius R g ~ 0.5d. For even smaller polymers 
the amplitude of the deformation becomes comparable 
to the polymer size. For polymers much larger than the 
interlamellar distances, our approach would need to be 
completed by accounting also for the pressure exerted by 
the polymer on the neighboring membranes (of order of 

k B Td- 3 ) g. 



Our treatment of the polymer induced deformations 
does not account for the fluctuations of the membrane. 
Thermal fluctuations arise spontaneously in membrane 
systems on the scales above the de Gennes-Taupin per- 
sistence length: £ = a'cxp(^Y) (42|. When the elastic 
constants are of the order of k B T, the lenght £ is in the 
range 10 — lOOnm, and fluctuations become important. 
At the level of the first order perturbation scheme im- 
plemented in this paper, the pressure is an external field 
that only couple linearly to the deformation of the mem- 
brane, so the fluctuation spectrum is not perturbed. It is 
therefore important to extend to second order this type 
of perturbative expansion in order to determine both the 
corrections to the elastic parameters of the membrane 
and the corrections to the fluctuation spectrum itself. 

Acknowledgments: We gratefully aknowledge in- 
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Obukov and T.A. Witten. 



V. DISCUSSION AND CONCLUSION. 



We have developed in this paper a new picture to de- 
scribe the interactions between a grafted polymer and a 
flexible membrane. 

We have shown that the polymer behaves in fact as a 
small pressure tool. The pressure applied by this tool is 
very high close to the grafting point and decays sharply 
over a distance of order R g , the gyration radius of the 
chain. Through scaling arguments and numerical sim- 
ulations we further confirmed that excluded volume in- 
teractions mainly change the range of the pressure field, 
while its amplitude and functional form are rather inde- 
pendent of solvent conditions. 

Under the pressure of the polymer, a flexible interface 
assumes a characteristic deformation: each flexible sur- 
face is a pressure sensor. We have shown that for fluid 
membranes, the deformation has a pinched conic form. 
The exact shape of the pinch depends on the boundary 
conditions imposed upon the membrane, but its form at 
the center of the deformation field is rather universal. 

For many grafted polymers, the deformation field in- 
duces an interaction potential between the polymers. We 
have shown that this potential is attractive for polymers 
grafted on the same side of a membrane and repulsive 
for polymers grafted on opposite sides. This interaction 
potential rises interesting possibilities. For instance, by 
changing the temperature one might expect to control the 
aggregation behavior of polymers grafted on the mem- 
brane. Also, if many polymers are added to a vesicle, the 
attraction might bring many polymers to the same site, 
increasingly catastrophically the pressure until a deco- 
rated micelle or small vesicule detaches from the surface. 
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FIGURES 



FIG. 1. The pressure applied by a grafted polymer to the 
surface. We chosed here a = 0.1R g and arbitrary units of 
pressure. The insert stresses the scaling form r~ 3 close to the 
grafting point. 

FIG. 2. Comparaison of the calculated 

pressure and monomer concentration at the wall extracted 
from Monte-Carlos simulations, a) Freely-hinged chain of 200 
monomers, b) Freely-hinged chain of 200 monomers with ex- 
cluded volume. The continuous line is the expression (^) of 
the pressure. Note that there is no adjustable parameter for 
the gaussian chain. For the chain with excuded volume, the 
size of the polymer is R 2 = 1.5N 2v a 2 , with v = 0.6. 

FIG. 3. The pinched form of a fluid membrane under the 
pressure of a grafted polymer. The deformation profile is the 
actual calculated form for a freely standing membrane. 

FIG. 4. The membrane-mediated interaction potential be- 
tween two gaussian polymers, grafted to the same side of a 
membrane. Polymers grafted on opposite sides have the same 
functional form with the reverse sign. 

FIG. 5. The two-body interaction potential with the es- 
timated inter-chain repulsion. The plots show 2 ^lp as a 
function of l/R g for different values of the bending rigidity 
(in ksT units), with fixed parameter b — l/(27r) . 



FIG. 6. The second virial coefficient as a function of 
the bending rigidity (in ksT units), with fixed paramater 
b = 1/(2-71-). Chains aggregation is induced by negative B, 
which occurs at k < 0.6fcsT for this particular value of b. 

FIG. 7. The conformation of a star anchored to a plane. 
Inside each "blob" of size £(r) ~ r/ _1,/2 the branches have a 
single-chain behaviour. 

FIG. 8. The deformation profile of a fluid membrane for 
different strengths of the harmonic potential. The curves are 
plotted for a = 0.5, 1 and 2, the largest deformation corre- 
sponding to the weaker potential. 

FIG. 9. Amplitude of the pinch at the grafting point. The 
parameter a controls the strength of the potential. For a 
soft potential, ho varies like the logarithm of a. The insert 
emphasizes the scaling form for large values of a. 
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